Tomonaga-Luttinger Liquid Renormalized by a Single Impurity 
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Effects of a single impurity potential on bulk properties of a spinful Tomonaga-Luttinger (TL) 
liquid is studied. A boundary bosonization technique is developed to include the impurity potential 
of arbitrary strength V . Our new bosonization formula for fermion field t/; smoothly connects the 
two existing expressions in the strong [V — oo) and the weak (V = 0) impurity limits. With use of 
the formula, we found the TL parameters determined from the long-distance correlation functions 
are renormalized due to the partial transmission through the impurity potential. 
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Interactions between particles introduce collective fea- 
tures into many-body systems. The dramatic many-body 
effects are reported in a variety of one-dimensional (ID) 
systems such as quantum wires Jj] , carbon nanotubes ^ 
neutral atomic gases of fermions Q and bosons 0, 
in ID optical lattices. A single-particle excitation loses 
its weight by exciting infinite numbers of particle-hole 
pairs. Low energy physics of interacting ID particles are 
described by the collective excitations of density fluctu- 
ations: Tomonaga-Luttinger (TL) liquid @, 0, 13- Only 
the collective features emerge in the lowest energy scales 
for homogeneous systems. 

Once inhomogeneity is introduced to the system by a 
localized potential V{x) ^ VS{x), such as an impurity, 
the single-particle excitations come into play due to their 
backscattering, and strongly affect the lowest collective 
excitations. They cause the zero-bias anomaly in the ID 
electron transport • Effects of a single impurity in a 

TL liquid have been intensively studied in the context of 
the scaling of impurity potential 11 , . Roughly speak- 
ing, for a spin independent system, the electron transmis- 
sion probability scales toward 1 or as lowering temper- 
ature, depending on whether the interparticle interaction 
is attractive or repulsive. This is a prediction from the 
perturbative renormalization group (RG) analysis near 
the two fixed point values, y = or cx). The scaling flows 
at arbitrary V have not been accessible with some excep- 
tions, the model for which the exact solution is available 
at special interaction parameters [l^, or for weak inter- 



action limits near Fermi- liquid fixed points IJ, [15| . Re- 
cently, the numerical assessment was done for intermedi- 
ate region by a path-integral Monte Carlo method [3]. 
However, their result shows in some part non-trivial be- 
haviors in the scaling flows of the zero-bias conductance. 
Thus, the effect of a single impurity of arbitrary strength 
is still unclear after the pioneering works of 15 years ago. 

Theoretical difficulties in treating arbitrary V are con- 
centrated on the point that both the long-wavelength 
and short-wavelength physics are of equal importance. 
This bears a resemblance to the Kondo effect 17| and 



the Fermi-edge singularity in X-ray absorption spec- 
if. The difference from them is in the point 



trum 

that the local degrees of freedom and the collective de- 
grees of freedom originate from the same conduction elec- 
trons. Hence, a unified theory describing the single- 
particle field and long-wavelength collective TL field is 
needed. The aim of this paper is to construct the 
bosonization formula which gives the relation between 
the single-particle phase shift 9f and the bulk excitations 
of the collective mode [13, 21, HI]- The power law de- 



pendencies of various correlation functions are calculated 
with use of the obtained expression. We show explicitly 
the bulk excitations suffers a significant change due to the 
presence of a single impurity potential. The zero bias 
conductance is also determined for a given phase shift. 
We set h = ks = 1, the inverse temperature /3 — l/T, 
and the system size L. 

The Hamiltonian under consideration consists of three 
parts H = Ho-\- Hy -\- iJint, 



Ho = J dx ^Ptix) {pl/2m - £f) Mx), (1) 
Hy = V [ dx S{x)^bl{x)Mx), (2) 



da; / dx' ncr{x)Ua.cr'{x ~ x')na-'{x'), (3) 



Hi, 



where ipa and rio- = t/'J.V'ct are the electron field and den- 
sity of spin (T = ±l=t,4. Hq, iJv, and Hmt represent, 
respectively, the kinetic energy, an impurity potential, 
and the two-body interactions. 

Before treating the spinful model, let's first derive a 
bosonization formula for a spinless model. We start with 
the one-body scattering problem 



-dl/2m + V6{x)) 4>{x) = ek4>{x). 



(4) 



The solutions of even (P = +1) and odd {P = —1) parity 
are: '/)p=i.fe(a;) = ^2/Lcos {\kx\ — 6'f) and (/)p=_i^A;(a;) — 
\/2/L sin (kx) . The momentum k is defined in a semi- 
infinite space fc > 0, and the scattering phase shift 
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= tan ^{mV/kp) is approximated at the Fermi level. 
The fermion field is given by ip{x) = J2p k>o 0P,A:(a^)cp,fc, 
where cp^k (cp j.) is annihilation (creation) operator of 
a state (P,fc). In treating the two-body interaction in 
bosonization technique, it is convenient to use the de- 
composition 



-ikpx 



(5) 



This procedure is brought forward easily for the homo- 
geneous case. In this case, however, one must keep in 
mind that these two chiral fields depend on each other. 
To fulfill the constraint, we introduce a fermion operator 
cp^fc, which is defined in a full momentum space (— cxd < 
k < cxd) and satisfies cp_fc = cp,_A;(= cp.fc). In addition, 
we keep the anti-commutation relation {cp j,, cp'^fc'} = 

Sp,p'Sk,k'- The chiral field tpr (t = ±) consists of two 
field 

i^rix) = e-'-'^^^-^''-^+,,r{x) + e-'^^/^^^iA^),{6) 
where each of them comes from the even or odd functions. 



i'P^rix) = — e'^ 



('-'^^"£p,rfc. (7) 



The slowly varying density Jt{x) = ijjl{x)ipr{x) is given 
by J^'P'(x) = ijjl^^{x)-ijjp',T{x) as 



Jrix) 



ire 



V 

Z^p p, 

-E 



e''?^-*^''p^p'(^)jP^P'(q). (8) 



The four density operators corresponds to the four kinds 
of particle-hole excitations characterized by the parities 
of a particle and a hole wavefunction (see Fig. [T]) . The 
scattering phase shift enters the above expression as 



'.p 



(a;) = Jp,_p- X P X (sgn(2;)6lF - 7r/2). (9) 



It is worth noting that the phase shift appears only in 
the density operators with odd parity PP' = — 1. It is 
naturally understood from the fact that in order to ex- 
cite the density wave of odd parity, a relative number of 
fermions between the right and the left of the impurity 
must change, i.e. the particles must tunnel through the 
impurity potential. It is convenient to express the four 
kinds of bosonic density operators in terms of pseudo spin 
defined by Ct{x) = ('0+i,r(a;), '0-i,r(2:))*- The pseudo 
spin density is defined by (i — 0, 1, 2, 3) 



(10) 



where 70 is a 2 x 2 unit matrix, and (71, 72, 73) are Pauli 
matrices in the pseudo spin space. The pseudo spin den- 
sities satisfy the commutation relation, 



(11) 




Jr 



Jr 



FIG. 1: Four types of the density operator J^'^ represent a 
pair excitation of a particle with the parity P = ±1 and a hole 
with the parity P' = ±1. In the each diagrams, energy levels 
of the even (odd) parity states are shown in the right (left) 
side of the vertical axis (energy) where the particles (holes) 
are represented by filled (open dotted) circles. 



We rewrite _ffo + Hy in bilinear form of them by linearis- 
ing the spectrum around the Fermi level as 



TIVp 
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J2 Y.Jr{im-<i)- (12) 



Next, we will consider interactions between particles. 
We assume that the relevant physics originating from 
the particle correlations can be captured with a (32-54)- 
model [i^l which takes into account only the forward 
scatterings and neglects the backward scattering gi . This 
approximation is valid quantitatively as far as the Fourier 
components of the interaction potential in Eq. ^ satisfy 
UaM'{q ^ 0) UaM'il ^ 2A;f) as in the case of long- 
range interaction. Then, the interaction Hamiltonian can 
be written 

ffint = y^_y da; [giJr{x)Jr{x) + g2Jr{x)J-r{x)] 

= [9AJ%q)J'A-q)+92J'r{q)J\{-q) 

+ €Y^^(Qp){g^jl{q)jl{-q)^g2Jl{q)j\{-q)) (13) 

+ COs2(0f) {g^Jl{q)Jl{-q) - g2Jl{q)JlA-q))\ ■ 

In order to get the second equation, we kept only the 
terms nonvanishing after the integration over x in the 
r.h.s. of the first equation. The nonvanishing terms sat- 
isfy the parity conservation (P1P2P3P4 = +1). Accord- 
ing to the commutation relation in Eq. (Ill|) . the bosonic 
ladder operators are introduced as 

J\{q) ^ \qL/2l,\^'^ (e{TqMq)+e{-Tq)b\{-q)) . (14) 

By introducing a Bogoliubov transformation 



h{q) 



sinh(0,) cosh(0O ) \l3\{-q) ' '^ ^^ 



the total Hamiltonian is diagonalized as iJ = 
Yjq,i '"i\q\{Pi{q)Pi{q) + \)- The rotation angles depend on 
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dimensionless coupling parameters ^2(4) = .92(4)/27rwF, 



tanh(20o) 
tanh(202) 



-.92 



54' 



tanh(2(?;>i) 



-g2 sin^(0F 



32 cos 







1+54 cos2(0p) ' 



1+54 sin^(6'F) 
tanh(203) = 0. (16) 



The four-component TL liquid is characterized by four 
TL parameters Ki = e^"^' {i = 0, 1, 2, 3). The parameter 
for the charge Kq is exactly the TL parameter in a ho- 
mogeneous system. For clarity, we shall call K{= Kq) a 
bare TL parameter. We now reach our main result, the 
boson representation of the fermion field 



tprix) 

Xt{x) 



L ^ ~ 



(17) 



g,P,P' 



q>0 * ^ 



iq 



rboiq) + rsgn(a;) sm{0F)bi{q) 



cos(0F)fo2(g))-li.' 



(18) 



To get this expression, we put the aforementioned con- 
straints cp_/c = cp^-k i.e. bi{q) = bi{—q). It is easily 
checked that this formula recovers the known expression 
in the weak and strong barrier limits by putting^ — 
or 7r/2. For example, in the strong barrier limit [24|, we 
have 

1¥ 



Xrix) = r J — (^l/_ft" cosgx + iT\^ sin qx 

q>0 ' ^ 



X P' {sgn{x)q) + h.c. 



(19) 



where /3'(±g) = iPo{q) ± Pi{q)) / {iV2) . With our for- 
mula, we clearly see the operators for a; > and a: < 
are separable for 9f = t^/^. In this case, the TL liquid is 
regarded as completely disconnected two pieces. 

The extension of these results to a spinful model is 
straightforward when the potential barrier is indepen- 
dent of spin Vi = Vi — V, since charge (c) and spin (s) 
variables are completely decoupled in the TL Hamilto- 
nian. It is simply done by adding real spin variables 
a = ±1 =t,i- The system is described as an eight- 
component TL liquid with the Hamiltonian {H)q = 

J2^=o-3,JJ^^3M{0J.M0]■M + l)■ Following the con- 
ventions [23|, we denote the dimensionless coupling con- 
stant for the charge (spin) gn'^'' — {gn±g:^)/2 for n = 2, 4. 
TL parameters are defined Kc{s),i = e^*^ where (j)c{s),i 
is given by putting —> (j)c(s),i and g„ -> g^*"' in 
Eq. (fTB|) . In the same way, the phase field XT,cr = 
{Xt,c + ^Xt,s)IV^ of the chiral fermion ii)r,a oc e'"^^^-" 



is obtained by putting x-! 



Xt,c(s) 



and bi 



c(s) A 



Eq. (|T8|) . In the case of the spin dependent potential bar- 
rier, a mixing of spin and charge degrees of freedom oc- 
curs through the difference in the scattering phase shifts. 
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FIG. 2: TL parameters K*^^^ determined by Ea. (l22p are plot- 
ted as function of transmission probability ^ with for several 
values of the interaction parameter with — 



In that case, however, one can solve the problem byintro- 
ducing a rotation in the spin and charge space H, 2^ . 

So far, we see the spinful TL liquid with a single im- 
purity is described as eight-component TL liquid. This 
indicates the bulk theory is subject to change in the 
presence of a localized scatterer. This change can be 
evaluated from the difference in the correlation functions 
in the bulk region \x\ 3> v-pP. Similarly, the renor- 
malization of TL parameters is also evaluated from the 
scalingdimensions of the long distance correlation func- 
tions [8|. The correlation functions of charge density 
wave Ocdw = 4'1 (ri^-T^cT, spin-density wave (sdw), sin- 
glet superconductivity (ss), and triplet superconductiv- 
ity (ts) are calculated as {O'' {x)0{x')) cx \x — x'\~'^ for 



\x I S> vp(3. The exponents are 



^cdv 



Kss = Ac-1 + + Kts = Ac-1 + A^.,-!, 

where, for j = c, s and — ±1, 



(20) 



2A, 



{K,,oY + sin^{ep){Kj,iY -f cos2(0p)(Xj- 2)- 



By comparing Eq. (pOI with those of homogeneous sys- 
tems the renormalized TL parameters {K*^ K*) of an 
effective two-component TL liquid must satisfy 



{K:r + {Kir - ^c,. + A,y 



(21) 



for arbitrary combinations oiu^v' = ±1. Generally, map- 
ping the eight-component TL liquids to two-component 
model is impossible. However, a solution of a set of two 
equations derived from Eq. pT|) . 



K* 




(22) 



gives the reasonable value since it is exact at least to the 
linear order in the small coupling constants g' s, 

K* = l-[(^^-^^\fj{9i+9i)- (23) 
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Thus, we accept K* in Eq. ([22)) as the effective TL pa- 
rameter. K* is plotted as a function of ^ in Fig. [2l The 
result indicates that the TL parameter, which character- 
izes the bulk quantity of the lowest energy scale in a TL 
liquid, is renormalized due to the presence of a single 
impurity. Since our calculation starts from a scatter- 
ing problem, all the single-particle states are under the 
effects of the single impurity pot ential and implicitly in- 
cludes the Friedel oscillation IJ] which has the singular 
effects on the electronic correlations especially in one di- 
mension. It is worthwhile to note that K* approaches 1 
for the intermediate value of ^ in Fig. [2l This indicates 
that the partial transmission (or the partial reflection) 
destabilizes the two-body interaction effectively than in 
the homogeneous system {\Kj^Q — l^jiO ~ 1|)- This 
could be understood as follows. The impurity scattering 
causes the phase mismatch between the density waves 
through Eq.®. Thus, there occurs a decrease in the 
overlap integral of the density waves due to the presence 
of the impurity scattering. This results in the weakening 
of the two-body interaction. While the TL parameters 
depend on the strength of the impurity scatterings, the 
zero temperature phase diagram, which tells the most 
dominant correlation functions, remain the same as in a 
clean system [23]; From Eq. ([20|l . the phase boundaries 
are determined by Ac.+i ~ ^c.-i and Ac,+i = ^c,-i, 
which is Kc,Q = ^s,o = 1- 

Finally, we refer to the transport properties. The 
charge current is defined by jc(a;, t) — dtXr.cix, t)/2Tr. 
The zero bias conductance G is calculated with the use 
of Kubo formula and Eq. (fT8|) . Since a DC electric field 
which is symmetric E{x) — E{~x) and concentrated 
around the impurity only couples with the current of a 
charge mode of i = 2, we obtain 



G = 




(5| + gi)cos2(0F) 



(24) 



This result shows resemblance to the Landauer formula 
for non- interacting quantum wire G/Gq — cos^(0f) with 
difference of the prefactor. However, it should be noted 
that the fermion phase shift given by Eq. (jH) is exact 
only for noninteracting systems, since the phase shift, 
which is determined by the impurity potential V, should 
be modified by the electronic correlations. RG equation 
of the phase shift may be evaluated from the correlation 
functions of boundary operators, such as the local den- 
sity of states [2^. One can reproduce the RG equation 
valid for arbitrary transmission but for the weak inter- 
Further discussions on the renormalization 



action 

of the transmission amplitude will be shown somewhere 
else. 

In conclusion, we have studied the bulk properties of 
a TL liquid with a single impurity potential using our 



new bosonization method. Our theory relates a one-body 
scattering phase shift to the properties of the collective 
excitations. The new bosonization formula smoothly con- 
nects the existing formulae valid for dp = or tt/2. TL 
parameters, which characterize the long distance corre- 
lation functions, are found to be renormalized due to the 
presence of a single impurity. 
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